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SUMMARY. 

Is  examined  the  problem  of  the  aerodynamic  rotation  of  the  plane 
of  the  motion  of  the  flight  vehicle,  which  realizes  deorbit  of  ISZ 
[MC3  -  artificial  earth  satellite]  into  the  dense  layers  of  the 
atmosphere,  with  its  subsequent  escape  to  assigned  flight  altitude 
above  certain  fixed/recorded  by  point  earth's  surface.  Are  indicated 
the  sufficiently  simple  methods  of  the  construction  of  the 
approximate  program  of  control  by  roll  attitude,  that  makes  it 
possible  to  satisfy  limitations  to  finite  values  of  phase 
coordinates . 
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Introduction. 

The  study  of  a  question  about  the  optimal  control  of  the  orbital 
apparatus,  which  accomplishes  three-dimensional/space  maneuver  in  the 
atmosphere,  is  of  interest  in  connection  with  the  fact  that  the  use 
of  lift-drag  ratio  on  the  apparatuses  of  this  type  makes  it  possible 
to  considerably  widen  the  possibilities  of  their  use/application, 
also,  in  a  number  of  cases  to  substantially  decrease  the 
general/common/ total  power  expenditures,  required  for  the  rotation  of 
the  plane  of  the  motion  of  apparatus,,  in  particular  on  the  rotation 
of  the  plane  of  initial  orbit  [1]. 

The  efficiency  of  the  use  of  aerodynamic  forces  for  the 
accompl isiiment  of  the  maneuver  of  orbital  apparatus,  until  now,  was 
studied,  as  a  rule,  in  connection  with  the  problem  about  the  rotation 
of  the  orbital  plane.  Power  expenditures,  required  for  the 
accomplishment  of  maneuver,  were  evaluated  according  to  the  value  of 
ideal  velocity,  necessary  for  the  compensation  for  aerodynamic  drag. 
It  was  established/installed,  that  the  aerodynamic  method  of  the 
rotation  of  the  orbital  plane  proves  to  be  energetically  more 
advantageous  in  comparison  with  rocket-dynamics  at  the  angles  of 
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rotation,  which  exceed  Ai=7°-9°,  and  as  the  approximate  law  of 
control  of  attitude  of  roll  7  with  a  sufficient  accuracy  can  be 
accepted  program  7=const. 

In  the  present  work  the  aerodynamic  maneuver  of  orbital 
apparatus  (Fig.  1)  in  the  initial  phase  of  trajectory  of  deorbit, 
which  is  characterized  by  the  presence  of  intense  reflections  from 
the  dense  layers  of  the  atmosphere,  is  investigated.  In  this  case  the 
use  of  engine  thrust  is  excluded. 
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Key:  (1).  Terrestrial  track  of  initial  orbit.  (2).  Terrestrial  track 
of  trajectory.  (3).  initial  orbit.  (4).  conditional  boundary  of  the 
atmosphere . 

Page  4 . 

FORMULATION  OF  THE  PROBLEM. 

In  spherical-high-speed  coordinate  system  without  taking  into 
account  the  rotation  of  the  Earth  the  equations  of  motion  of 


S 
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apparatus,  which  possesses  lift-drag  ratio,  take  the  following  form 


where 


■^  =  V'sine; 
at 


"  _  .■  cosiigvl 

dt  Vcosa  L  2  r  J 


do 

l/cos9slmj 

W 

r 

dk 

V'C0S®C0Si| 

dt 

”  r  cos? 

c.S 

_ 

5  — .  »  P 

Po« 


-ZA 


(la) 


(1 


r  =  /?  +  A.  /?- 6371,21- 10* 

ft  =  9,81  M.cek*,  z  —  ^; 

Tj  -  heading/course  angle  between  the  projection  of  velocity 
vector  on  the  plane  of  the  local  horizon  and  the  parallel. 


(p,  X  -  geographic  latitude  and  longitude/length. 


7  -  attitude  of  roll  of  the  apparatus  (positive  value  7 
contributes  to  an  increase  of  77). 


Let  us  note  that  the  motion  of  apparatus  in  the  plane  of 
scanning/sweep,  described  by  subsystem  (la),  does  not  depend  on  the 
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angle  of  bank. 

Let  us  assume  to  the  piecewise-cont inuous  control  ^(t)  the 
limitation  of  the  form 

Tmin  Taiii.  ) 

is  superimposed,  and  the  functional,  maximized  on  the  final 
diversity,  isolated  in  the  phase  space  by  equations 

X(f,)-X>=0;  (3) 

IS  certain  value 

9 (A),  -niti))-  (4) 

The  control  7{t),  which  maximizes  1(7)  with  the  the  free  t^,  and 
also  the  free  moment/torque  t,  of  the  descent  of  apparatus  from  the 
orbit,  we  will  call  optimum,  and  the  corresponding  phase  trajectory 
x(t),  which  satisfies  conditions  ( 3 ), -extremal . 

Page  5. 

APPROXIMATE  PROGRAM  OF  CONTROL  OF  ROLL  ATTITUDE. 

We  will  use  conventional  assumptions  [2],  connected  with  the 
assumption  of  the  smallness  of  angle  Bit)  and  the  smallness  of  flight 

altitude  h(t)  in  comparison  with  a  radius  R  Earth; 

sin  6^6;  co8  8~1;  r~/?. 


(5) 
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Passing  to  the  new  independent  variable/alternating/variable 
angular  flying  range  s,  with  the  help  of  relationship/ratio 
ds=(V/R)dt  let  us  convert  system  of  equations  of  motion  (1)  taking 
into  account  (5)  to  the  following  form: 

V~-a,pV,  j 


9  =  a,  p  cos  7  — 


(6a) 


h  —  RH-, 

t|  =  a,  ?  sin  f  —  cos  T)  tg  ?; 
^  —  sin  T|; 

i  —  SSUL- 

cos?  ’ 


where 


= 


■  (4r=*)- 


Entering  in  accordance  with  the  principle  of  L.  S.  Pontriagin's 
maximum  [3],  [4],  let  us  register  the  function  of  Hamilton  H  and  the 
conjugated/combined  system  of  equations: 


Xj,  7)-(p»cosT  +  ;7,sinT)a,p— /^v,a,p  V+/7*/?6-f 


itafaL 
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'2gR 

Pv-Pv‘^x9  —  P*-^> 

P*  =  —PhR> 

?*==■  —Pv  <^x  Vi?  +  (p%  cos 'T  4-  />,  sin  t)  ?r , 

,  s'nij 

Pi  =  — /7,sln7]tga  -/’»cosT„ 


;>,=;»,  COST) 


1 


cos-  9 
px  =*0. 


-Pi 


COS  T)  sin  ? 
cos*  9 


(^) 


Optimal  control  7  on  open  set  of  the  allowed  values  is 
determined  from  conditions  [3]: 


s«nx 


ojt 


where 

Page  6 . 


—  A*  ~P* 

— p====  ;  COS-Iaf»  —  ^ — J-— — j-  » 

^  Pi  +  P*  ^P»  +  Pl 

. .  pi  COS  T.pt  -f  Pi  sta  Ton  <0.  ^ 


(9) 


Taking  into  account  that  in  the  problem  H(Pirxi»Tj™d,  R=const 
question,  from  (7)  and  (8)  it  is  not  difficult  to  obtain  the 
relationships/ratios  following  equivalent  between  themselves: 

.  -  .  P»-}-Pp»  =  Ct-.  V  -V  '  (^0)'  • 


where 


pj=-^p*/?z® ij  29.  • 

g  —  —  PiCOSi^9  4- pv  -f-  Ptsinn,  Rz  —  const. 


(12)  . 
(13) 


in 
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Equation  (11)  in  the  range  of  the  velocities,  which  do  not 
exceed  circular  (v<v'gR) ,  describes  the  process  of  forced 
oscillations. 

In  connection  with  this  function  (-s).  and  consequently 
[according  to  (10)],  and  f unct ior  they  oscillate  about  the  zero 

values.  Therefore,  using  certain  averaged  value  of  f unct ion  p»(s),  it 
is  possible  to  obtain  the  approximate  program  of  control  of  roll 
attitude . 


Actually/really,  if,  for  example,  value  P*  is  determined 
condition  (12),  assuming/setting  in  the  first  approximation, 
constant. 


P%ep' 


const. 


—  1 


(14) 


f  rom 
PKiS) 


then  formula  for  the  approximate  program  of  control  of  roll  attitude 
will  take  the  following  form: 

_ El _  (15) 

Comparison  with  the  formula 

of  that  being  determining  the  optimum  program  of  roll  attitude  in  the 
mode/conditions  of  quas i-stat ionary  gliding/planning  (6=0)  shows  that 
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the  presence  in  denominator  of  expression  (15)  of  value  gives 

further  information  about  a  change  in  the  instantaneous  value  of  the 
flight  path  angle  d{s).  This,  in  turn,  gives  the  further  possibility 
to  actively  affect  the  flight  trajectory  due  to  the  appropriate 
assignment  in  (15)  of  the  /alue  of  free  parameter 

It  is  characteristic  that,  in  contrast  to  (15),  formula  (16) 
defines  attitude  of  roll  7  at  those  points  in  the  trajectory,  in 
which  the  second  derivative  of  value  considered  as  the  function  of 
argument  x=ln  V/V, ,  takes  zero  values. 


Page  7. 


Actually/really,  function  pt{x)  with  v<V^  satisfies  with 
respect  to  the  variational  equation  of  the_form 


dI  4.  fiK  t  ^ 

ai,"  (v-  -37  • 


(17) 


and,  therefore,  at  points  in  the  trajectory,  at  which /7('=0,  value 


is  determined  by  the  relationship/ratio 

Use  in  formula  (15)  of  the  first  integrals  of  system  (6) 
obtained  by  V.  F.  Illarionov  (8)  in  the  form 
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Pf  =  ~  +  pi  cos  (X  -  >./). 

Pi  —  — ^-0  “  Pf  sin  (^  —  ^)  cos  f  ^  pi  sin  ? 


(lb) 


makes  it  possible  to  express  the  instantaneous  value  of  the 
approximate  roll  attitude  in  the  function  of  phase  coordinates  and 
some  constants  of  integration  pi,  pi.  p(.  p*cp,  determined  from  the 
boundary  conditions  of  problem  at  the  j- end/lead  of  the  trajectory. 
The  need  for  the  integration  of  interconnected  circuit  in  this  case 
completely  is  excluded. 


One  should  say  that  during  the  determination  of  the  approximate 
law  of  control  of  attitude  of  roll  according  to  formula  (15),  valid 
for  open  set  of  the  allowed  values  7,  limitation  (2)  can  be  performed 
via  usual  "truncation"  of  value  7. 

FURTHER  SIMPLIFICATIONS.  STRUCTURE  OF  THE  APPROXIMATE  CONTROL  OF  THE 
MODULUS  OF  ROLL  ATTITUDE  IN  THE  CLASS  OF  PIECEWISE  CONSTANT 
FUNCTIONS. 

In  spite  of  sufficient  simplicity  of  formula  (15),  the 
establishment  of  characteristic  properties  in  the  structure  of 
optimal  control  of  roll  attitude  nevertheless  remains  difficult.  In 
connection  with  this  let  us  pass  to  the  examination  of  the 


I 
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approximate  control  in  the  class  of  the  piecewise  constant  functions. 

Let  us  consider  the  case  of  the  single  insert  ion/ immers ion  of 
apparatus  in  the  atmosphere  at  first.  We  will  count  the  initial  orbit 
of  equatorial.  Then  latitude  o  in  (1)  will  characterize  the 
instantaneous  value  of  the  value  of  parallax  relative  to  initial 
orbit,  the  execution  of  the  second  of  conditions  (3)  can  be  provided 
with  the  appropriate  selection  of  the  moment/torque  of  the  descent  of 
apparatus  from  orbit  t,- 

Will  introduce  into  the  examination  the  conditional  boundary  of 
the  atmosphere,  which  corresponds  to  certain  height/altitude 

(19) 

higher  than  which  motion  can  be  considered  Keplerian.  Let  us  take 
condition  (19)  as  the  sign/criterion,  determining  end  point  of  time 
— 4  Of'  the  ascent  path,. 


Using  formulas  of  Keplerian  theory,  the  not  difficult  initial 
problem  of  reducing  to  the  problem  of  maximization  on  the  boundary  of 
the  atmosphere  (at  moment/torque  t=t^)  of  the  functional  of  the  form 


when  the  limitation  of  the  form  is  present. 
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Thus,  let  us  consider  the  atmospheric  trajectory  phase,  in  which 
we  will  assume/set  the  effect  of  aerodynamic  forces  predominating 
over  all  remaining  forces,  which  function  on  the  apparatus,  and  we 
will  consider  that  the  angles  and  ©  of  three-dimensional/space  turn 
in  the  limits  of  one  immersion  into  the  atmosphere  are  small: 

cosij^l;  sio?~f;  cos9~l'.  (22) 

Under  these  conditions  the  system  of  equations  (6)  takes  the 
following  form: 


8  =  a,  .0  cos  7,  i 

(23a) 

h  -  /?0;  ' 

T  =  a,psinT,  j 

b 

?  =  1 

(23(y 

>:=i.  1 

We  will  assume  that  the  scalar  function  ©(V(t,),  0(ti),  t?  ( t , )  ) 
in  (4),  determined  at  end  point  of  time  ti,  is  such,  that  the 
maximization  of  the  final  velocity  V(t,)  in  the  constant/ invar iable 
ones  S(ti)  and  i}(ti)  continuously  corresponds  to  the  maximization  of 
functional  1(7).  In  that  case  as  this  is  not  difficult  to  show,  the 
solution  of  the  initial  problem  about  the  maximization  of  functional 
(20)  under  the  condition  of  connection/communication  (21)  is  found 
among  the  solution  set,  which  correspond  to  the  problem  about  the 
maximization  of  the  velocity  at  the  end  of  atmospheric  section  I/,  at 
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height/altitude  A  =  /i*  with  different  fixed  values  of  final  parameters 
*»•  >  and  ?*•  In  connection  with  this  the  unknown  solution  can  be  found 
with  the  method  of  the  parametric  optimization  of  values  0^.  n*  and  ?». 

We  will  use  this  fact  and  let  us  consider  the  case  of  "small" 
optimum  (required)  angles  AI* (t;,,  (pj  of  the  aerodynamic  rotation  of  the 
plane  of  motion,  which  do  not  exceed  in  the  modulus/module  of  the 
corresponding  angles  of  rotation  Ai»,  obtained  in  the  problem  about 
maximization  at  height/altitude  =  in  free  parameters  and  ?*• 
constant/invariable  parameter  and  constant  attitude  of  roll  7. 

Thus , 

(24) 

In  this  case  reaching/achievement  of  fixed/recorded  parameters 
■»1*  and  ?*  can  be,  obviously,  provided  without  the  further  expenditures 
of  energy  due  to  the  appropriate  programming  of  the  angle  of  the  bank 
(exchange  of  sign  7  does  not  affect  motion  in  the  plane  of 
scanning/sweep) . 

Taking  into  account  of  the  aforesaid,  value  and  in  (6)  it  is 
possible  to  unfasten  and  subsequently  to  examine  only  system  of 
equat ions 
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9  =  a,  p*  cos  7,  j  (25) 

P  =  —  Rz^p.  1 

Let  us  note  that  since  the  trajectory  of  apparatus,  which  is 
characterized  by  reflection  from  the  dense  layers  of  atmosphere 
(d>0) ,  interests  us,  we  will  assume  that  the  value  maximum  on  the 
modulus/module  of  roll  attitude  does  not  exceed  certain  limiting 
value  lTnp«i!‘ 


(26) 
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Accepting  as  the  argument  the  monotonically  increasing  variable 

6(s),  we  convert  system  (25)  to  the  form: 

d\n  V  _  1 

d6  A  cos  7  ’ 

df  _  ze  '  (-"J 

dH  a,  A  cos  7  ‘  ; 

As  a  result  we  come  to  the  following  task;  it  is  necessary  to 
determine  the  control  y  on  set  of  allowed  values  (26),  which 
maximizes  functional  Ini/*  at  the  fixed  final  values 
variable/alternating/variable  ji?  and  equal  to 

P.-P*.  8.  =  9*. 

Initial  parameters  V,,  0,  and  p,  are  assigned. 

According  to  the  principle  of  L.  S.  Pontriagin's  maximum,  the 
function  of  Hamilton  and  the  conjugated/combined  system  of  equations 
take  the  following  form: 

H -sr) 

A  =  0;  =  (29) 

Constant  in  (29)  value  Pj,  which  corresponds  to  coordinate  In  V, 
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let  us  assign  as  follows: 

p,  ss  const  <  0. 


(30) 


Then  optimality  condition  of  control  7,  which  consists  of  the 
minimization  of  Hamilton's  function^  (7),  takes  the  form: 


llLia  ^  (®)>0- 


(31) 


Key:  ( 1 ) .  with. 


where 

(32) 

there  is  a  function  of  the  changeover  of  control  7. 


To  the  moment/torque  of  the  changeover  of  boundary  I7I 
corresponds  certain  value  of  argument,  determined  from  the  condition 
of  equality  to  zero  functions  of  changeover  of  the  formula 


Pi  2 


(33) 


From  formula  (33)  taking  into  account  (30)  it  follows  that 


r.) 

|<0  p,<0, 

\>0  npH 


(34) 


Key:  (1) .  with. 
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If  we  now  take  into  account  that  to  the  condition  of  the 
positiveness  of  the  function  of  changeover  it  corresponds: 


l<Cr^npH 


(35) 


Key:  (1) .  with. 


that  of  (30),  (34)  and  (35)  easily  is  established  that  independent  of 
Ipil  the  minimum  on  the  m.odul  us /module  roll  attitude  is  optimum  on 
the  ascent  path  [fl(s)ii0]  when  Pj<0.  But  if  Pi>0,  then  attitude  of 
roll  txi-lTf.u,  is  optimum  respectively  on  the  descending  branch  of 
trajectory  [^(s).<:0]. 

We  will  consider  that  lTLia  =  0.  Then  taking  into  account  the 
uniqueness  of  the  moment/torque  of  the  changeover  of  limiting  angle 
of  bank  from  the  aforesaid  easily  is  establ ished/ instal led  the 
following  structure  of  optimal  in  (25)  control  of  the  modulus/module 
of  attitude  of  roll  (Fig.  2): 


ill 


Pt<0. 

jhl«« 

.  1  t09«  1  —  )  g 

npH  9,<9(sX9,; 
8,  <9(s)<9»; 

9.<0, 

(36> 

Pt>0, 

9X9(s)<9,; 

9*<9(5X9j. 

K>o, 

(37) 

Key :  ( 1 ) .  with . 


In  each  specific  case  the  optimality  of  one  or  the  other  type  of 
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control  is  determined  by  the  boundary  conditions  of  task. 


It  is  possible  to  show  that  the  physical  sense  of  parameter  pj 
the  task  in  question  is  in  connection  with  the  following: 


Pi  =  — 


^^rr.ix  (Pi) 


(38) 


In  view  of  the  linearity  of  system  of  equations  (27)  function 
/m»x(?»)  is  convex.  Hence,  for  example,  it  follows  that  in  the  region 
of  the  solvability  of  task  (Fig.  3),  determined  (at  the  fixed  value 
of  angle  by  the  inequality 

?*  mio  <?»<?*  mil*  (39) 

in  S6Qrn6nt.  Pit  mat  I  there  is  a  unique  point  ?»  =  ?»•■  at  which 

functional  on  set  of  permissible  controls  U  reaches  absolute 

maximum  [(/„„(?*,)— sop  / (f)]. 


1/4  i 


iiV 


V:- 

1  i 

1  i 

Fig.  2, 


Fig.  3. 


Fig.  2.  Structure  of  optimal  roll  control  during  single 
insertion/immersion  of  apparatus  in  the  atmosphere. 

Fig.  3.  Region  of  solvability  of  task. 

Page  11. 


Therefore 


C" 

<^P*  |<0  npH 


Key ;  ( 1 ) .  with. 


Condition  (40)  taking  into  account  (38)  determines  the 
dependence  of  the  sign  of  value  p,  and,  consequently,  also  the 
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structure  of  optimal  control  7  from  the  boundary  conditions  of  task. 


For  determining  the  limiting  values  of  final  density  P»  mjt  and 
Pkmio  at  fixed  value  9»,  it  is  necessary  respectively  to  maximize  and  to 
minimize  in  (27)  coordinate  ?»  at  the  free  final  flight  speed.  With 
the  help  of  the  principle  of  L.  S.  Pontriagin’s  maximum  we  obtain  the 
appropriate  programs  of  the  control: 


P*  P*  mini 


P*  P*frii» 


(§> 

unpti  I  8(s)<0. 

0  8  (s)>-  0; 


(P 

npH  8(s)<:0, 
mpH  6  (s)>  0. 


(41) 


(42) 


Key :  ( 1 ) .  with . 


Value  P*«  corresponds  to  the  value  of  final  density  P*-  obtained 
as  a  result  of  the  solution  of  the  problem  about  the  maximum  of  value 
In  in  free  parameter  P*  and  fixed/recorded  8^.  in  this  case  we  have 
p,=0,  and  therefore  when  P*=P»» 

iToptl^O.  (4.3) 

Integrating  system  (27)  under  the  laws  of  control  (41)-(43), 
taking  into  account  (38)  and  (40)  we  come  to  the  following 
inequalities,  which  isolate  the  regions  of  the  boundary  conditions  of 
tasks,  which  correspond  to  sign-constancy  of  parameter  p^  (Fig.  4): 


1 
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— —  1  -C- 


p* 


a4, 


(4S> 


Pi>0. 


Fig .  4  .  Fig .  5 . 


Fig.  4  . 
?  2  .  1  - 
control 


To  determination  of  regions  constancy  of 
program  of  the  control  of  fox'm  (41),  2  - 
of  form  (43),  3  -  program  of  the  control 


«*COi  T 


tftsp  — 


apc« 


€  ak  ’ 


«,s«rctf  -r - 


BpCA 


sign  of  parameter 
program  of  the 
of  form  ( 42 ) , 


Fig.  5.  Dependence  oiC^)  for  two  types  of  optimal  control - optimal 

control  with  p,<0  form  (36),  -  the  same,  with  Pj>0  form  (37), 


s  »retg 
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By  the  use  of  inequalities  (44)  and  (45)  is  determined  the  sign 
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of  parameter  pj  for  any  values  ?»  and  9^,  ass igned  in  the  region  of  the 
solvability  of  task,  and  therefore  the  corresponding  type  of  the 
optimum  program  of  control  of  the  modulus/module  of  attitude  of  roll 
[ { 36)  or  ( 37)  ] . 

It  is  characteristic  that  when  ?*  =  ?,  in  each  of  the  regions  of 
boundary  conditions  (44)  and  (45)  are  valid  the  relat ionsh ips/ rat ios 
of  form  and  19*1<1®«I  respectively.  In  connection  with  this  the 

program  of  control  of  the  attitude  of  roll  of  form  (36)  (Fig.  5) 
corresponds  to  the  more  intense  flight  of  apparatus  from  the  dense 
layers  of  the  atmosphere  (to  high  values  of  final  angle  l®»ii 

Let  us  note  that  on  the  basis  of  the  physical  sense  of  parameter 
Pj,  to  condition  p2<0  corresponds  decrease  with  decrease?*  when 

•n  — const,  or,  that  nevertheless,  decrease  with  increase^  when 

— const  Hence  the  sense  of  mutual  (converted)  task  lies  in  the  fact 
that  in  fixed/recorded  parameters  V*  and  P*  value  9*  with  pj<0(pj>0) 
should  be  maximized  (to  minimize).  It  is  analogous  in  fixed/recorded 
parameters  V»  and  9»  with  pj<0  (pj>0)  should  be  minimized  (to  maximize) 
value  P»-  The  latter  is  shown  in  Fig.  3  by  arrows/pointers . 

Fig.  6  depicts  the  typical  dependence  of  the  finite  values  of 
the  parameters  of  trajectory,  which  correspond  to  flight  altitude 
h‘=100  km,  and  also  apogee  altitudes  of  trajectory  9,  from  the  flight 
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path  angle  determined  at  the  height  taken  as  the  conditional 
boundary  of  the  atmosphere,  equal  to  A»  =  70  k.u.  The  value  of 

lift-drag  ratio  is  accepted  equal  to  K=1.5,  and  (sign 

7=const ) . 
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Fig.  6. 

Fig.  6.  Effect  *»  on  parameters  of 
flight  altitude  h=h‘  under  law  of 

Key:  (1).  m/s. 


Fig.  7. 

trajectory,  which  correspond  to 
control  of  form  (36). 


Fig.  7.  Structure  of  optimal  roll  control  during  double 
insert  ion/ immers ion  of  apparatus  in  the  atmosphere. 

Key;  (1).  Outer-atmosphere  trajectory  phase. 
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Page  13. 

It  is  characteristic  that  as  a  result  of  an  increase  in  the 
duration  of  flight  at  the  maximum  roll  attitude  vith  increase/growth 
the  angle  of  rotation  of  the  plane  of  motion  Ai,  and  also  the  value 
of  lateral  deviation  o  continuously  grow. 

Using  the  results  obtained  above,  let  us  construct  the 
approximate  program  of  the  piecewise  constant  control  of  the 
modulus/module  of  roll  attitude  for  the  case,  when 

(46) 

For  this  purpose  in  the  examination  of  the  complete  system  of 
equations  of  motion  (6)  we  will  use  the  specific  character  of  the 
programs  of  control,  close  to  the  optimum  ones,  in  the  following  two 
tasks,  which  are  special  cases  of  the  task  in  question: 

a)  the  task  about  reaching/achievement  of  fixed/recorded  angle 
with  the  maximum  final  rate  in  free  parameters  ’i*  and  (free  angle 
Ai*); 


b)  the  task  about  the  rotation  of  the  plane  of  motion  to  preset 
angle  Ai*  ("n*- ?»)  with  the  maximum  final  rate  in  free  parameter  6,. 
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Specifically,  taking  into  account  that  in  the  second  task  as 
this  was  shown  V.  P.  Plokhikh,  the  program  of  the  attitude  of  roll  of 
form  7=const  on  the  functional  is  sufficiently  close  one  to  the 
optimum,  let  us  represent  the  general  routine  of  control  in  the  form 
(36)  and  (37),  assuming  that  the  minimum  value  of  modulus/modul e  I7I 
in  (2)  corresponds  to  certain  value  iVl  /  necessary  for  the  rotation 
of  velocity  vector  to  angle  =  i Ax»optl  —  i  ‘• 

Thus,  as  a  result  of  replacing  limitation  (2)  by  the  inequality 

(47) 

the  approximate  program  of  roll  attitude  in  the  class  of  the 
piecewise  constant  functions  takes  the  following  form: 


Pi  C  0, 

lTl  =  j 

hi 

e,<9(s)<e,; 

8,  <8  (S)  <  9»; 

».<0; 

(48) 

P»>0, 

1,1  =  1 

hi  ^’"npH 
hl.uMnpH 

8,<8(s)<6,; 

8,<9(s)<9». 

8,>0; 

(49) 

Key;  (1) 

with. 

In 

connection 

with 

task  "b" 

in  (43)  we 

have  9, —9^ 

and  in 

correspondingly,  9,  =  8,. 

One  should  say  that  since  the  value  of  angle  usually  is 

previously  unknown,  when  SAj,  ^0  the  parameter  I7I  in  (47)  is  subject 
to  optimization. 
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Let  us  note  finally  that  from  the  linearity  of  the  system  of 
equations  of  motion  (27)  the  analogy  in  the  structure  of  control  of 
the  modulus/module  of  roll  attitude  easily  follows  in  the  atmospheric 
trajectory  phases  in  the  case  of  the  double  i  r.sert  ion/ immers  ion  of 
apparatus  in  the  atmosphere.  In  this  case  the  type  of  the  approximate 
program  of  control  [(36)  or  (37)],  adjusted  as  a  result  of  the 
examination  of  the  task  about  the  maximization  of  the  flight  speed  at 
the  end  of  each  of  the  atmospheric  trajectory  phases  at  the  fixed 
values  of  the  angles  of  arrival  ®*iopt  and  depends  on  the 

value  of  the  latter.  Thus,  Fig.  7  shows  the  structure  of  control  in 
connection  with  the  case,  when  |A‘»op«t<lA<»l. 

In  the  examination  of  the  real  motion,  described  by  a  precise 
system  of  equations  (1),  it  is  convenient  to  use  the  fact  that  the 
program  7=^  is  dividing  two  types  of  the  approximate  in  (23)  control, 
which  is  characterized  either  by  the  amplification  of  the 
fluctuations  of  flight  altitude  [formula  (36)],  or  by  damping  the 
fluctuations  [formula  ( 37 ) indicated] . 

Page  14. 


In  this  case  as  the  fundamental  characteristic  of  the  approximate 


DOC  =  83087402 


PAGE  >5^ 


control  7  apogee  altitude  of  the  outer-atmosphere  ellipse  of 
trajectory  ^  =  (0*.  V"*,  p,),  can  be  accepted  since  when  p*=-const,  (in 

view  of  freedom  of  movement  in  the  plane  of  scanning/sweep  from  the 
angle  of  bank)  parameters  8^^  and  t'*  are  wholly  determined  by  the 
structure  of  control  of  the  modulus/module  of  roll  attitude.  Then, 
for  example,  in  the  case  of  the  single  insert  ion/ immers ion  of 
apparatus  in  the  atmosphere  when  should  be  used  a 

structure  of  the  control  of  form  (36),  and  when  — 

respectively  form  (37).  It  is  characteristic  that  if  in  the  first 
case  to  limitation  of  the  type  of  equality  h(ti)=h‘  to  final  flight 
altitude  h(ti)  is  equivalent  limitation  of  the  type  of  the 
conditional  inequality  of  form  h(ti)?h‘,  then  in  the  second  case  the 
limitation  of  form  h(ti)<h‘  is  equivalent  to  it. 


EXAMPLE 

Let  US  consider  the  task  about  the  rotation  of  the  plane  of 
motion  with  the  yield  condition  at  the  end  of  the  maneuver  into  the 
point  (X*,  <p‘,  hM,  with  the  maximum  final  rate. 

The  initial  parameters  of  motion  we  will  consider  the  following 


h(0)=130  km;  V(0)=7820  m/s;  d{  0)=-l'’,5;  ")(0)=0. 
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As  the  limitation  to  the  roll  attitude  let  us  take  inequality  of 
l7l$87“. 

For  the  construction  of  the  program  of  control  in  the  class  of 
the  piecewise-cont inuous  functions  we  will  use  formula  (15). 

Taking  into  account  (^i)  =  0  with  the  designations  of  the  form 

I  ccstsln? 
cos; cos 3 

1  sin ;  P)tcf  ^ 

from  relat ionsh ips/ rat ios  (13),  (15)  and  (18)  we  obtain 

I _  ^  cos;  sin  3  —  5  cos -cos  ^ 

D  cos ;  sin  ?  —  £  cos ;  cos  ?  -f-  sin  ;  ’ 

where  =  cos  ?  sin  (X  - 1'). 

fl  =  tg  s' cos  ?  cos  (X  — 1')  —  sin  ?.  .  ,7. 

£>  =  sin  ?  cos  T,  sin  (X  —  x')  -f  sin  n' cos  (X  -  x“), 
f  —  cos ? cos T]  -f  tg?'  (sin? COST] cos  —  X‘)  —  sin tj sin  (X  —  X')j. 

As  a  result  the  program  of  roll  control  proves  to  be 
biparametr ic .  The  parameters  ^  and  ^  are  selected  from  the  condition 
of  falling  in  point  (o',  hM. 

In  the  case  in  question  is  realized  the  double 
insertion/immersion  of  apparatus  in  the  atmosphere  (Fig.  8),  and 
change  7  leads  to  the  amplification  of  the  fluctuations  of  flight 
altitude  h,  since 


I _ 1 _ 

ly'p;  {t,)+pUt,)~(PHc,Rr 
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Aj  I  ^  0 .  A.=o<CA’.  (51) 

Let  us  compare  the  now  obtained  program  of  control  with  the 
approximate  control  7  (t)  in  the  class  of  the  piecewise  constant 
functions.  Taking  into  account  (51),  let  us  consider  the  program  of 
the  control  of  form  (48),  assuming/setting  attitude  of  roll  7  by 
identical  for  each  of  the  sections  of  the  insertion/immersion  of 
apparatus  in  the  atmosphere.  As  a  result  we  come  ,to  the  four-stage 
approximation  of  function  7(t):  iTil  =  ITj  1  =  It  lm«»  Ta  “=  7*  =  T-  The  unknown 
parameters  are  modulus/module  |Vl  »  and  also  two  moments/torques  of 
changeover  It  Lax  to  value  7=^« 

Page  15. 


From  Fig.  8  it  is  evident  that  the  obtained  law  7(t)  on  the 
whole  correctly  reflects  the  behavior  of  stepless  control.  The  error 
in  the  value  of  the  maximized  final  rate,  obtained  as  a  result  of  a 
similar  approximation,  composes  the  value  of  order  1%.  It  is  obvious 
that  the  accuracy  of  approximation  can  be  increased,  if  we,  for 
example,  7,  and  7*  optimize  independent  of  each  other. 

It  should  be  noted  that  in  spite  of  a  larger  number  of  free 
parameters,  the  construction  of  the  piecewise  constant  program  of 
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control  indicated  proves  to  be  simpler  in  comparison  vith  the 
procedure  of  the  use  of  formula  (50).. 
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